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ABSTRACT 

We report observations of three rotational transitions of molecular oxygen 
(O2) in emission from the H2 Peak 1 position of vibrationally excited molecular 
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hydrogen in Orion. We observed the 487 GHz, 774 GHz, and 1121 GHz lines 
using HIFI on the Herschel Space Observatory, having velocities of 11 km to 
12 km and widths of 3 km s~^. The beam-averaged column density is N(02) 
= 6.5x10^^ cm^^, and assuming that the source has an equal beam filling factor 
for all transitions (beam widths 44, 28, and 19"), the relative line intensities imply 
a kinetic temperature between 65 K and 120 K. The fractional abundance of O2 
relative to H2 is 0.3 - 7.3x10"^. The unusual velocity suggests an association 
with a ~ 5" diameter source, denoted Peak A, the Western Clump, or MF4. 
The mass of this source is ~ 10 Mq and the dust temperature is > 150 K. Our 
preferred explanation of the enhanced O2 abundance is that dust grains in this 
region are sufficiently warm (T > 100 K) to desorb water ice and thus keep 
a significant fraction of elemental oxygen in the gas phase, with a significant 
fraction as O2. For this small source, the line ratios require a temperature > 180 
K. The inferred O2 column density ~ 5x10^^ cm~^ can be produced in Peak A, 
having N(H2) — 4x10^^ cm^^. An alternative mechanism is a low- velocity (10 
to 15 km s~^) C-shock, which can produce N(02) up to 10^^ cm~^. 

Subject headings: ISM: molecules; molecules - individual (oxygen); Orion molec- 
ular cloud 



1. INTRODUCTION 

1.1. Chemistry of Oxygen in Dense Interstellar Clouds 

Oxygen is the third most abundant element in the universe, and its form in different 
regions is thus of considerable astronomical importance. Whittct (2010) has recently reex- 
amined the question of oxygen abundance as a function of density in the interstellar medium 
(ISM), and finds that a large fraction of oxygen is unaccounted for at the higher densities 
which generally correspond to the molecular ISM. Models of gas-phase chemistry that do not 
treat the condensation of species as icy mantles on grains have predicted molecular oxygen 
(O2) to be abundant in well-shielded regions, with X(02) = N(02)/N(H2) greater than 10~^, 
and thus the most abundant molecular species after H2 and CO, with predicted abundance 
exceeding that for water (Herbst & Klemperer 1973; Prasad & Huntress 1980; Graedel, 



Herschel is an ESA space observatory with science instruments provided by European-led Principal 
Investigator consortia and with important participation from NASA. 



-4 - 



Langer, & Prerking 1982; Leung, Herbst, & Huebner 1984; Langer & Graedel 1989; Bergin, 
Langer, & Goldsmith 1995; Marcchal et al. 1997b). The conditions that generally (but not 
in every model) maximize the abundance of O2 are low metal abundance (consistent with 
depletion on grains), high densities, and sufficient time for oxygen chemistry to come into 
steady state. 

The gas-phase chemistry of O2 in dense clouds is relatively simple. Formation of O2 
is initiated by the reaction of atomic oxygen with Hj^ to form 0H+. }if is produced by 
cosmic ray ionization of H2 and subsequent reaction of the resulting with H2. The OH"*" 
reacts with H2 to produce H2O+, and by a second reaction, HsO^. This ion dissociatively 
recombines with an electron, and the two channels of interest result in H2O and OH. The 
latter species can react with an oxygen atom to form O2. 

Uncertainties in several of the key reaction rates have been reduced by laboratory mea- 
surements. The dissociative recombination of HsO^ has been measured using heavy-ion 
storage rings by Vejby-Christensen et al. (1997), by Jensen et al. (2000) and by Neau et al. 
(2000). The reaction rate for all three of these measurements is in reasonable agreement in 
magnitude and temperature dependence, with a large thermal rate coefficient a ~ 3xl0~^ 
cm^s~^ at 50 K, and a temperature dependence oc T"^-^. These results are consistent with 
the earlier results by Mul et al. (1983) obtained using a different technique. The branching 
ratio has been a subject of some uncertainty, the papers cited above all put the fraction of 
the recombinations that result in production of OH in the range 0.5 < /(OH) < 0.7. Thus, 
the effective rate coefficient for formation of OH by dissociative recombination of HaO"*" at 
the temperature of dense molecular clouds is about 2x10^^ cm'^s^^. 

The rate at low temperatures of the reaction OH + O — )■ O2, the final step in the produc- 
tion of molecular oxygen, has been the subject of some controversy. Theoretical calculations 
by Harding et al. (2000) were in moderately good agreement with the measurements for 
temperatures between 150 K and 500 K (Howard & Smith 1981; Smith & Stewart 1994), 
and indicated the peak value of the rate coefficient to be A; = 5.5x10"^^ cm^ s~^ at ~ 80 
K, with a factor ~ 2 reduction at 10 K and 1000 K. Measurements by Carty et al. (2006) 
sampled a range of temperatures, but these authors recommended use of a single value, k — 
3.5 ± 1.0x10-^^ cm^ s"^ over the temperature range 39 K < T < 142 K. 

The calculations by Xu et al. (2007) suggested a rate that while only somewhat lower 
than that of Carty et al. (2006) at temperatures above 80 K, fell dramatically at lower 
temperatures. This opened the possibility that the O2 abundance in cooler clouds would 
be much lower than predicted by the above astrochemical models. However, more elaborate 
quantum mechanical calculations by Qucmener, Balakrishnan, & Kendrick (2008) showed 
a relatively weak temperature dependence, and a rate coefficient at 10 K that is 70 times 
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larger than predicted by Xu et al. (2007). 

Recent quantum calculations by Lique et al. (2009a) give rate coefficients that are a 
factor ~ 2 greater than the measurements of Carty et al. (2006) over their temperature 
range. The rate coefficients from these latest theoretical calculations do fall at temperatures 
below 30 K, but even at 10 K are > 3xl0~^^ cm^ s~^, thus guaranteeing that the OH + 
O — > O2 + H reaction will not be a low-temperature bottleneck for O2 production. The 
agreement between measured and theoretical calculations of the reaction rates suggests that 
the O2 production rate is reasonably well understood and confirms the purely gas-phase 
chemistry prediction of large X(02) in all but possibly the coldest clouds. 

The expectation of a high abundance of O2 had previously been an impetus to include 
this species in calculations of gas-phase cooling (Goldsmith & Langer 1978), and it was 
found to be moderately important at intermediate densities, despite its weak magnetic dipole 
transitions. The major obstacle to studies of astronomical O2 is the strong absorption by this 
same species in the Earth's atmosphere. The large scale height of O2 compared to that of H2O 
means that even going to high altitude sites does not allow for observations of Galactic O2 in 
the interstellar medium, where only the energy levels producing millimeter and submillimeter 
rotational transitions will likely be populated. Although absorption observations of UV 
transitions are possible from space (Smith et al. 1984), clouds with sufficient column density 
to have detectable transitions of O2 are likely to have sufficiently large dust opacities that 
such observations will prove difficult. 



1.2. Grain Surface Oxygen 

A possibility that must be considered is that O2 is resident on dust grain surfaces. Solid 
O2 has been searched for in interstellar clouds through its weak feature at 6.45 /xm, with 
negative results (Vandenbussche et al. 1999). The most stringent abundance limits come 
from analysis of the solid ^^CO line profiles in low- and high-mass star-forming regions, which 
indicate solid O2 concentrations of less than 10% with respect to water ice (Boogert, Blake, 
& Tielens 2002; Pontoppidan et al. 2003). Taking a typical water ice abundance of 10^^ 
relative to H2, this indicates sohd 02/gas-phase H2 abundance ratios < 10"^ within a factor 
of 2. Thus, at most 10% of elemental oxygen is locked up in O2 on grain surfaces. 

Laboratory measurements show that solid O2 ice is only slightly less volatile than CO 
ice with a binding energy of ~ 900 K, and thermally evaporates at 16-18 K under typical 
interstellar conditions (Acharyya et al. 2007). If trapped in water ice, a small fraction of 
O2 can desorb at higher temperatures together with H2 (CoUings et al. 2004). Once on 
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the grains, solid O2 is chemically active. A number of recent laboratory experiments have 
demonstrated that O2 is rapidly hydrogenated by atomic H to form OH, H2O2 and H2O 
in ices, even at temperatures as low as 12 K (Miyauchi et al. 2008; loppolo et al. 2008, 
2010). Some gas-phase O2 can also be produced by photodesorption of H2O ice at high 
temperatures (Oberg et al. 2009). The nondetection of solid O2 on grain surfaces is thus 
reasonable, and overall, the amount of oxygen on dust grains is not a major contributor to 
the overall oxygen budget in dense interstellar clouds. 



1.3. Searches for Moleculcir Oxygen 

Several different techniques have been utilized in the search for interstellar O2. In 
general, these normalize the derived column density (or upper limit thereto) to that of 
standard tracer, generally an isotopologue of carbon monoxide, from which the molecular 
hydrogen column density can be derived. The resultant fractional abundance of or upper 
limit to the fractional abundance of O2 is obtained assuming that the distributions of the 
species are the same. The isotopologue ^^O^^O has levels and transitions that are not allowed 
for the homonuclear common isotopologue, and some of these are at frequencies at which the 
atmosphere is reasonably transparent, notably the (N,J) = (2,1) — j- (0,1) transition at 233.95 
GHz (Black & Smith 1984; Marechal et al. 1997c). Ground-based observational searches for 
16QI8Q j-^avg been carried out with nondetections giving upper limits X(02) somewhat lower 
than those for CO (Goldsmith et al. 1985; Liszt & Vanden Bout 1985; Liseau et al. 2010). 
A ground-based search (Pagani, Langer, & Castets 1993) yielded a tentative detection in 
one source that was not confirmed by further observations (Marechal et al. 1997a) with 
significantly improved limits. Fuente et al. (1993) searched for the 234 GHz transition of 
^^O^^O in three dark clouds, and measured la upper hmits N(02)/N(C0) < 0.15, 0.19, and 
0.29 for in TMC2, L134N, and B335, respectively. 

The common isotopologue of molecular oxygen can be observed using ground-based 
telescopes from sources having velocities sufficient to Doppler-shift the observed frequency 
"out from under" the pressure-broadened terrestrial O2 absorption. The lower rotational 
transitions of O2 referred to here are shown in Fig. 1. This technique has been used to 
search for O2 in the galaxy NGC 7674 by Liszt & Vanden Bout (1985), using the redshifted 
119 GHz line, who found X(02) < 1.4x10"^ (1 a). Goldsmith & Young (1989) searched for 
the 119 GHz O2 transition in the galaxy Vll Zw 31, determining a 1 a upper limit X(02) < 
4x10"^ Combes et al. (1991) searched for O2 in NGC 6240 as weU as ^^O^^O in two Galactic 
sources, obtaining upper limits X(02) < 1x10"^. Liszt (1992) added three galaxies having 
recession velocities between 8800 and 12,000 km and carried out additional observations 
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of NGC 7674, finding N(02)/N(C0) < 0.05 at the la level. Combes & Wiklind (1995) 
searched for red-shifted 368 GHz and 424 GHz O2 lines in the BL Lac object B0218+357 
having z = 0.685, determining N(02)/N(C0) < 0.014. Combes, Wiklind, & Nakai (1997) 
observed this same source, but in the 56.3 and 119 GHz O2 transitions, and obtained an 
improved upper limit N(02)/N(C0) < 2xl0~^. Frayer et al. (1998) searched for O2 in 
Markarian 109, obtaining an upper limit N(02)/N(C0) < 0.31 in this low-metallicity galaxy. 
All of the above ground-based observations had angular resolutions < 1'. 

Observations of O2 are possible from balloon altitude for sources within the Milky Way 
with sufficiently large velocities. Olofsson ct al. (1998) searched for the 424 GHz transition 
of O2 from an altitude of 39.5 km with a 60 cm telescope giving a beam width ~ 5'. They 
did not detect O2, but were able to determine upper limits in two sources, NGC7538 and 
W51, with N(02)/N(C0) < 0.1. Since X(CO) is expected to be ~ lO"^ in warm regions, 
these hmits correspond to X(02) < 10~^. 

The best environment for searching for molecular oxygen in the Milky Way as well 
as other galaxies is clearly above the Earth's atmosphere, from a spacecraft. Out of a 
number of proposals, two modest-sized space missions dedicated to submillimeter wavelength 
spectroscopy have been successfully developed and operated. The first of these was NASA's 

Submillimeter Wave Astronomy Satellite (SWAS) launched in 1998, which observed the 487 
GHz O2 line (Melnick ct al. 2000) with a beam width 3'. 5 x 5'. Numerous clouds in the Milky 
Way were observed, but there were no definitive detections of O2. Upper limits reported by 
Goldsmith et al. (2000) were as low as X(02) < 2.6x10"^, and a combination of data from 
9 sources yielded (X(02)) = [0.33±1.6(3ff)] x 10"^. 

The Odin Satellite (Nordh et al. 2003) was launched in 2001. It was equipped with 
a larger antenna and more fiexible receiver system than those of SWAS. Observations of 
O2 were primarily of the 119 GHz transition with a beam width of 9' (Pagani et al. 2003), 
giving upper limits to X(02) in half of the sources being as low as 10^^. This survey included 
relatively warm GMCs, but due to the low energy of the upper level of this transition, the 
upper limits for O2 in dark clouds were significantly lower than those obtained with SWAS. A 
study of the Sgr A cloud in the Galactic Center yielded a Scr upper limit X(02) < 1.2x10"^ 
(Sandqvist et al. 2008). Wilson et al. (2005) used Odin to search for O2 in the Small 
Magellanic Cloud galaxy, and after 39 hours of integration obtained a 3a upper limit X(02) 
< 1.3x10^® in this source that has a (total) oxygen abundance less than 0.2 of that in the 
Orion nebula in the Milky Way. 

A very deep SWAS integration on a mosaic of positions in the p Ophiuchi cloud yielded 
a detection of the 487 GHz line at modest statistical significance (Goldsmith et al. 2002). 
The 5.5 km s~^ velocity of the observed emission line fell within the red-shifted wing of a 
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Fig. 1. — Lower rotational energy levels and transitions of O2. Each level is described by 
its rotational quantum number N and total angular momentum quantum number J, with J 
= N-1, N, N+1. This designation applies rigorously only to Hund's case (b), as described 
by Lique (2010), but is generally used to designate O2 levels. The frequencies of the 
allowed transitions are given in GHz, together with the spontaneous decay coefficients in 
s~^. Transitions with AN = are indicated by dotted lines; these have frequencies in the 
vicinity of 60 GHz with the exception of the N,J = 1,1 to 1,0 transition at a frequency of 
118.8 GHz. The seven transitions that fall in the frequency range covered by Herschel HIFI 
are shown in color, with the three relatively strong transitions discussed here (487, 774, and 
1121 GHz) are shown in heavy dark blue. See Marechal et al. (1997a) for additional details. 
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low- velocity molecular outflow. Since the column density of the outflowing gas is modest, the 
fractional abundance of O2 derived was quite large, X(02) = lxlO~^. This result was not 
confirmed by subsequent Odin observations of the 119 GHz line, which did detect an emission 
feature, but at a velocity of 3.5 km s~^, and linewidth of 1.5 km s~^, corresponding to the 
characteristics of the quiescent cloud rather than the outflow. The line is quite weak and due 
to the large column density of the bulk of the source, the derived fractional abundance of 
O2, X(02) = 5xl0~^ (Larsson et al. 2007), is relatively low. The consensus emerging from 
these observations of O2 is that the abundance of this species in a variety of different types 
of molecular clouds is dramatically lower than predictions of standard gas-phase chemistry. 
Instead of an abundance X(02) ^ 10~^, the upper limits are as low as 10~^, only slightly 
above the level of the Odin p Ophiuchi detection. 

However, even preceding the arrival of data from space missions, there had been theo- 
retical models of the ISM that predicted lower O2 abundances. Xie, Allen, & Langer (1995) 
investigated turbulent diffusion, which in principle could accelerate mixing of material from 
UV-irradiated cloud edges with the material deeper inside. The result is reduction in the 
abundances of certain molecules including O2. A model employing circulation has also been 
developed by Chieze & Pineau des Forets (1989), which can result in a significant reduction 
in the abundance of O2 if the characteristic circulation time is as short as 10^ yr. Le Bourlot, 
Pineau des Forets, & Roueff (1995) studied a bistabihty in interstellar chemistry that under 
certain conditions can result in a high ionization phase in which the abundance of O2 is re- 
duced by a factor ~ 100 relative to standard gas-phase chemistry. This is largely restricted 
to regions of low-to-moderate density, but can be favored by high cosmic ray fluxes. 

Results from SWAS and Odin indicated that the low abundance of O2 is widespread, 
and included very well-shielded regions of high density with n(H2) > 10^ cm~^. Additionally, 
the abundance of H2O was found to be substantially below that expected from gas-phase 
chemistry models. To address this, Bergin et al. (2000) developed a model invoking depletion 
of atomic oxygen onto dust grains, and conversion to water ice remaining on the grain 
surface. This model results in an increase of the gas-phase C/0 ratio, which has the effect 
of dramatically reducing the fractional abundance of O2 deep inside clouds. Hollenbach et 
al. (2009) later showed that at intermediate depths {Ay ~ 2 — 8), photodesorption of water 
ice from grains maintained a sufficiently high gas-phase abundance of elemental oxygen that 
abundances X(02) ~ 3xl0~^ were achieved in these regions, independent of the gas density 
or incident ultraviolet flux. 

The situation regarding O2 has remained sufficiently puzzling that an Open Time Key 
Program (Herschel Oxygen Project, or HOP) was proposed and selected to carry out a 
survey of warm molecular clouds in three rotational transitions of this elusive species. This 



- 10 - 



paper presents the results for some of the observations in Orion, which we give in §2. In §3, 
we discuss the identification as O2. In §4 we discuss collisional excitation of O2, the source 
temperature and O2 cohimn density derived assuming the source fills all HIFI beams. In §5 
we analyze the effects of beam filling and source offset, and discuss the column density and 
fractional abundance of O2 that apply to different possible source geometries. We discuss 
shocks and warmed grains as sources of the observed O2 in §6. In §7 we present a summary 
of our results and conclusions. 



2. HIFI Observations 
2.1. Data Taking 

These data were taken with the Heterodyne Instrument for the Far Infrared (HIFI) 
(de Graauw et al. 2010) on the Herschel Space Observatory (Pilbratt et al. 2010). The 
observations were carried out on operational days (OD) 291, 297, 474, 504, and 505 in 
March, August, September, and October 2010. The position observed is the direction of 
the NW maximum of H2 ro-vibrational emission (Beckwith et al. 1978). The coordinates 
arc aj2ooo = 5'*35™'13'^.7, 5j2ooo = -05°22'09". Double beam switching was employed. For 
observations of Orion the spacecraft pointing constraints meant that the reference positions 
were located within 10° of an east-west line, with one position being 3' east of the source 
position and the other being symmetrically located to the west. Eight settings of the local 
oscillator were used to enable sideband deconvolution. HIFI band la was used for the 487 
GHz observations, band 2b for the 774 GHz observations, and band 5a for the 1121 GHz 
observations. Observations of the J = 7 - 6 C^^O transition were obtained simultaneously 
with the 774 GHz O2 data, and the results are discussed in §3.2. 

The total (on-source + off-source) integration time for the 487 GHz observations was 
440 minutes for H polarization and 398 minutes for the V polarization, due to the failure 
of one segment of the observations for the latter. The integration times were 430 and 409 
minutes for the 774 GHz and 1121 GHz observations, respectively. The FWHM beam width 
is 44" at 487 GHz, 28" at 774 GHz, and 19" at 1121 GHz, while the main beam efficiency, 
eMB, is equal to 0.76, 0.75, and 0.64 at the three frequencies, respectively (Olberg 2010). 
These long integrations provided a test of the performance of the HIFI radiometers. As 
discussed in Appendix 1, the noise in the final spectra is a factor ~ 1.6 greater than expected 
theoretically. Figure 2 shows the HIFI beams on an image of the H2 emission in Orion, and 
indicates some of the compact sources in the vicinity. 
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Right Ascension (J2000) 



Fig. 2. — Image of vibrationally excited 2 fim H2 emission in Orion (from Bally et al. 2011). 
The coordinates of H2 Peak 1 towards which Herschel was pointed for these observations are 
aj2ooo = 5'^35'"13^7, 5j2ooo = -5°22'09". The three circles indicate the FWHM beam sizes 
at the three observed frequencies: 44" at 497 GHz, 28" at 774 GHz, and 19" at 1121 GHz. 
The position of the Hot Core (aj2ooo = 5^35""14^5, 5j2ooo = -5°22'33.6") is indicated by 
the black circle 10" in diameter, approximately the size of this region (Wilson et al. 2000). 
The Peak A - Western Clump - MF4 - Cnt D source (see references in text) is indicated by 
the yellow circle, 5" in diameter. The infrared source IRc2 (Downes et al. 1981) is indicated 
by the black cross. 
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2.2. Data Processing 

The data were processed using the standard HIFI pipehne in HIPE 4.10 (Ott 2010), 
then exported into FITS format for further analysis using the IRAM Gildas software package 
(http://www.iram.fr/IRAMFR/GILDAS). The HIPE task FitHifiFringe was used as part of 
the pipeline processing to remove standing waves from the data. This was followed by the 
task DoDeconvolution which utihzed the multiple local oscillator settings to extract single 
sideband spectra. In Figure 3 we show data for a hmited range of velocities around the 
frequency expected for each of the lines. The continuum as well as the wings of nearby 
strong spectral lines have been removed by baseline fitting. A broader view of these extremely 
sensitive submillimeter spectra will be presented elsewhere. 

The two orthogonal hnear polarizations of each of the HIFI bands are received by sepa- 
rate feedhorns and processed independently. They thus inevitably have slight relative point- 
ing offset on the sky, which due to the complex nature of the Orion region, gives us useful 
information on the source location, which is discussed in §3.3. 

We see that there is a feature for each transition that peaks between 11 and 12 km 
s~^(all velocities refer to the Galactic Local Standard of Rest). The results from fitting a 
single Gaussian to each of these features are given in Table 1. The presence of numerous 
other lines, especially close by to the 487 GHz O2 line means that the uncertainties in the 
line parameters are increased by the uncertainties in the fitted baselines. For the lowest- 
frequency transition, this may well be the dominant contribution to the total error, so that 
the uncertainties given in Table 1 must be regarded as conservative. The 487 GHz transition 
shows a second component near 5 km s~^, which is a separate entry in Table 1. We attribute 
this feature primarily to O2 emission from the Hot Core being picked up only in the large 
beam at this relatively low frequency. This is the only transition in which there is evidence 
for this second velocity component. The issue of the O2 signal coming from the Hot Core, 
as well as other possible contributions, is discussed further below. 

3. Identification as Molecular Oxygen 

Given that the integrated intensities represent 8.5 o", 22 cr, and 5 a statistical features, 
there is little doubt that we have detected emission lines at the frequencies of each of the 
three O2 transitions that we have observed at the H2 Peak 1 position. The question of 
whether the features are due to molecular oxygen depends on the agreement and plausibility 
of the velocities, the consistency of the line intensity ratios, and the presence or absence of 
other carriers at the three frequencies. 
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Fig. 3. — Spectra of the three rotational transitions of O2 observed with the Herschel HIFI 
instrument towards the Orion H2 Peak 1 position. The intensity scale is the main beam 
brightness temperature (antenna temperature divided by main beam efficiency). The red 
solid lines are fitted Gaussians, which all peak at LSR velocities between 11 and 12 km s^^. 
For the 487 GHz transition there is clearly a second feature, indicated by the broken red 
line, which we attribute primarily to O2 emission from the Hot Core at Visr = 5.5 km s~^, 
as discussed in the text. 



Table 1. Parameters of Observed O2 Lines 



Transition 


Frequency 


/ Tmbdv 


Vo 


AV (FWHM) 


Tmb (peak) 




(GHz) 


(K km s-i) 


(km s^-*^) 


(km s~^) 


(K) 


3,3-1,2 


487.249 


0.095 ±0.011 


10.96 ± 0.16 


3.05 ±0.40 


0.029 






0.113 ±0.013 


5.70 ± 0.23 


4.43 ±0.65 


0.024 


5,4-3,4 


773.840 


0.177 ±0.008 


10.96 ± 0.07 


2.91 ±0.16 


0.057 


7,6-5,6 


1120.715 


0.091 ±0.018 


11.87 ± 0.29 


2.87 ±0.62 


0.030 
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The ~ 11 km s"^ component of the 487 GHz transition and the two higher frequency 
transitions arc all consistent in terms of line width. The central velocity of the 1121 GHz 
line is higher by about 0.9 km s^^ than the central velocities of the other two lines. The 
formal uncertainty on the fitted central velocity is 0.3 km s^^, and the difference in velocities 
is thus three times this, although less than 1/3 of the FWHM line width. It is thus possible 
that, assuming all three lines are from O2, that the velocities are in agreement but we are 
seeing a statistical fluctuation in the fitting of the central velocity, most fikely of the highest 
frequency line. Alternatively the highest frequency line is sampling different gas. 



3.1. Interloping Lines 

It seems improbable that a single molecular species other than O2 would have three 
spectral lines that match so closely with those from the species of interest. It is more plausible 
that the three spectral features in question are produced by three different "interloper" 
molecules, but in this case they additionally must have essentially identical line widths. 
While we cannot rule out this possibility, it should be born in mind that in this region of 
Orion, there is a wide variety of line shapes and widths due to the presence of multiple 
kinematic structures, so the species responsible would have all to share one specific set of 
line profiles and widths, further reducing the likelihood of this explanation. 

We have carried out a targeted analysis to assess hkehhood that the observed spectral 
features are due to known molecular species. To do this, we first selected all molecular tran- 
sitions found in the SPLATALOGUE catalog (http://www.splatalogue.net) falling within 
5 km s~^ of any of the three O2 lines observed, and having upper state energies less than 
1000 K. There were two species (CH3NH2 and CH3OCHO) having lines possibly interfering 
with the 487 GHz O2 line, three species ((CH3)2CO, g-CH3CH20H, and OS^^O) close to the 
774 GHz line, and 6 species (CH3I3CH2CN, g-CH3CH20H, HCCCHO, CH3NH2, g-CsH, and 
t-CHaCHCHO) proximate to the 1121 GHz fine. 

For each species, wc used the CASSIS program (http://cassis.cesr.fr) to simulate the 
emission spectrum for a 5" diameter source having line width between 2 km s"^ and 4 km s^^, 
assuming LTE with kinetic temperatures between 75 K and 150 K. In every case except one 
(discussed below) we found that we had sufficient frequency coverage to include undetected 
stronger lines located somewhere in the HIFI spectrum. From the nondetections, we could 
set an upper limit to the column density of the species in question that was so low that 
all lines in the vicinity of the O2 transitions were significantly weaker than the observed O2 
main beam brightness temperatures. The range of frequency covered in the HIFI spectra 
ensured that transitions from states of comparable energy above the ground state to those 
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producing possibly interfering lines were observed. Thus, the gas temperature could not be 
adjusted in a manner to make the lines close to O2 relatively more intense than transitions 
of other carriers that were not detected. An example of this is (CH3)2CO (acetone). This 
molecule has a transition only 12 MHz from the 773.8397 GHz O2 line, but also has a group 
of transitions near 774.4 GHz predicted to be an order of magnitude more intense than the 
line close to O2. We do not detect the higher frequency lines. The sensitivity of our data 
is such that the possible contribution to the O2 line is more than an order of magnitude 
smaller than the observed feature. In this manner, we were able to eliminate 11 out of the 
total of 12 potentially interfering molecular species from consideration. 

One problematic species is GH3OCHO (methyl formate), whose 406 34 ^ 396 33 transition 
in the v = 1 state of the E species has rest frequency 487251.99996 MHz, compared to 
487249.3755 MHz for the 3,1 - 3,2 O2 transition. This is the only O2 transition with possible 
interference from this species. The upper level energy is 705 K. A number of narrow features 
in the full HIFI spectrum are reproduced assuming that methyl formate is present emitting 
with a centroid velocity of 8 km s~^ and a line width of 2.5 km s^^. If we assume that we 
are directly pointed at a 5" source having kinetic temperature 100 K, we derive a column 
density of 1.5x 10^^ cm^^. The exact value of the column density depends on the location and 
characteristics of the source relative to the HIFI beams, as well as the complex characteristics 
of emission from this species found by Favre et al. (2011). If the emission is entirely from 
the MF4 source (see Figure 2), the column density required would be considerably greater 
than found by Favre et al. (2011) so there may well be an extended emission component 
contributing. With these values defining the methyl formate emission, a significant portion 
of the 5.5 km s^^ feature of the 487 GHz O2 line seen in Figure 3 could be due to this species. 
Note that to reproduce the 11 km s^^ O2 feature, the CH3OCHO would have to be at 12.6 
km s^^, a value that is outside the range found by Favre et al. (2011), seen in our Herschel 
data, or found for any of the other species in this region of Orion. 

The 5.5 km s^^ feature seen in the 487 GHz O2 transition in the bottom panel of Figure 
3 is at the velocity observed for most emission from the Hot Gore region near KL. If all or 
part of this feature is due to O2, the fact that we see it only in the 487 GHz spectrum is 
not likely to be a result of excitation; as discussed at length in §4.1, the O2 transitions are 
likely to be close to thermalization and the relative intensity of the 774 GHz line compared 
to that of the 487 GHz line will not be very different from unity. We are more likely seeing 
a result of the very different beam sizes and the offset between the Hot Gore and H2 Peak 
1. The angular offset is 27" and the 1/e beam radii are 26.3", 16.6", and 11.4", at 487, 774, 
and 1121 GHz, respectively, yielding coupling factors to a point source 0.34, 0.07, and 0.003, 
respectively. The finite size of the Hot core will modify these results slightly, but it is still the 
case that when pointed at H2 Peak 1, the 487 GHz beam will couple moderately well to the 
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Hot Core, but the two higher frequency beams will have very small coupling, explaining why 
the 5.5 km feature is seen only in the lowest-frequency O2 line. This issue is discussed 
further in §5. 

Another aspect of the confirmation that the three lines shown in Figure 3 are indeed 
produced by O2 is the "chemical poverty" of the H2 Peak 1 position relative to that of KL 
and the Hot Core. This is illustrated by comparing the two spectra in Figure 4, which shows 
a dramatic reduction in the strength of almost all lines at H2 Peak 1 compared to KL. Many 
species drop by more than an order of magnitude in integrated intensity. The line crowding 
in Orion KL is one of the reasons we chose to avoid this source in our initital observations 
of O2 with HIFI. A small number of features are stronger, or at least more readily visible 
at H2 Peak 1. At this point, we have not been able to make secure identification of most of 
these. The conclusion is that finding an "interloping species" mimicking the three putative 
O2 lines in Figure 3 is definitely much less likely at H2 Peak 1 than it would be at Orion KL. 



3.2. O2 Line Velocities and Source Identification 

Studies of the central portion of the Orion Molecular Cloud in different molecular tracers 
have resulted in a number of compact regions being associated with the velocity range 10 - 
12 km s~^, despite this being quite distinct from the nominal velocity of the Hot Core (5 - 
6 km s~^) and of the Compact Ridge including KL itself (8 - 10 km s~^). In what follows, 
we briefly discuss these sources, with our goal being to show that associating the observed 
velocity range of the O2 lines with one of these regions is not unreasonable. 

Pauls et al. (1983) identified a source which they denoted "Peak A" in NH3 (3,3) 
emission, at location aBi95o — 5'*32"*46*.7, 5bi95o -05°24'24". Although the centroid velocity 
of this source is 2 km s~^, the centroid of the spectrum at offset position (-2", 2"), thus 
displaced toward our H2 Peak 1 position, is ^ 10 km s~^. Wilson et al. (2000) discuss 
a 10 - 12 km s^^ component of the (4,4) transition of NH3. These authors provide little 
information about its location, but it is apparently the northwest portion of the Hot core 
and could be an extension of Peak A. As such it is likely within 5" of the Hot Core nominal 
position, offset by ~ 27" from our H2 Peak 1 position. Plambeck & Wright (1987) see a 
relatively distinct peak in HDO in a velocity bin spanning 10.5 to 15.5 km s~^ at a position 
coincident with Peak A, which is 8" to the northwest of the Hot Core, and thus about 21" 
from H2 Peak 1. Masson & Mundy (1988) present a marginally spectrally distinct feature 
in the 12 to 14 km s""*^ range in HC3N J = 12 - 11, in the overall spectrum of the Hot Core, 
and find that emission at around 13 km s^^ peaks in the "Western Clump" which is ~ 5" 
west of IRc2, and essentially coincident with Peak A. However, the velocity range 9.5 km 



-17- 



15.0 1^ 1 , , , , ^ , ^ 1 ^ ^ ^ 1 ^ , : 1 : r 



^ 10.0 - 




250.8 850.4 850.6 850.8 851,0 

Frequency (GHz) 



Fig. 4. — Spectra taken at the CSO in January, 2011, while observing the NO molecule. The 
upper spectrum is of KL, while the lower spectrum is of H2 Peak 1. The components of the J 
= 5/2toJ = 3/2 transition of NO are indicated by cyan lines, transitions of CH3OH by red 
lines, of C2H5CN by green lines, of SO2 and ^^S02 by magenta lines, of NS by black lines, 
and of CH3OCHO by black lines. Identifications of transitions in the signal sideband are 
indicated by solid lines, and of transitions in the image sideband by dashed lines. The source 
velocity has been assumed to be 9 km s~^, so that transitions from the hot core, typically 
having velocities 5 km s~^ to 6 km s~^ are slightly shifted from their identifying lines. 
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s ^ to 11.5 km s ^ is also characteristic of the northern ridge cloud, observed in a number 
of species by Wright, Plambeck, & Wilner (1996). 

Murata et al. (1991) carried out a study of CS J = 2 - 1 and J = 3 - 2, and found a 
number of interesting structures, including numerous clumps which are quite prominent in 
the velocity range 9.9-11.4 km s~^ in their interferometric map of the higher frequency line. 

Most of these lie along a northeast - southwest line defining the compact ridge, but there 
are three additional (unnumbered) peaks. The first is just 2" southeast of H2 Peak 1, the 
second is about 16" south of our pointing position, and the third is 7.5" east and 20.6" south 
of our H2 Peak 1 position. This last peak is coincident with the Western Clump - Peak A 
feature. While no details are given about these peaks, they seem to be part of an arc-shaped 
ring of emission that lies just outside the blue lobe of the CO bipolar outflow mapped by 
Erickson et al. (1982). Murata et al. (1992) reported 3 mm continuum mapping of the 
central portion of Orion, in which a number of condensations were identified. Of greatest 
relevance to the present study are the "Cnt A" peak which coincides with the second CS 
peak of Murata et al. (1991), and "Cnt D", which is coincident with the third CS peak, and 
thus with Western Clump - Peak A. 

A recent high-angular resolution study of methyl formate (CH3OCHO) by Favre et 

al. (2011) indicates that at least two small clumps have a second kinematic component 
at velocities higher than the 7 km s^^ to 8 km s~-^ characteristic of most of the emission 
in the center of the Orion molecular cloud surrounding the Hot Core and the Compact 
Ridge. One of these clumps, denoted MF4, has its second velocity component, in the three 
different methyl formate lines observed, at between 10.2 and 11.4 km s^^, thus quite similar 
to the velocities of the O2 hues. This condensation agrees in position with the Western 
Clump - Peak A feature, confirming that this velocity range does seem to be characteristic 
of this compact source. The issue of the identification of a specific condensation with the O2 
emission will be discussed further in what follows and in §6. 

One line observed with Herschel HIFI simultaneously with the 774 GHz O2 hne is the 
J = 7 - 6 transition of C^^O. A hint of a peak at velocity greater than 10 km s~^ can be 
seen in the ^^CO J = 7 - 6 line observed by Pardo, Cernicharo, & PhiUips (2005), but in 
the vicinity of KL the profiles are dominated by the broad wings of the outflow. In Figure 
5 we show the spectrum of the line of the rare isotopologue observed with Herschel HIFI. 
The asymmetry in the narrow component of the emission is immediately apparent, and we 
also show a four-Gaussian component fit, which is the minimum required to reproduce the 
observed line profile reasonably well. We give the parameters of the components in Table 
2. The presence of the wide component and the asymmetric narrow component leaves the 
uniqueness of the very narrow narrow component uncertain. It is clear, however, that there 
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must be a component of the gas having central velocity in excess of 10 km s~^ at the H2 
Peak 1 position. In §4.3 wc use the C^^O emission to determine the column density of H2 
and the fractional abundance of O2. 

Another relevant molecule is NO, which can be produced by reaction of atomic nitrogen 
with OH, analogous to the pathway for producing O2. We have observed the complete set of 
the fine- and hyperfine-structure components of the ^Ili/2 J = 5/2 3/2 transition, which 
fall between 249.6 and 251.6 GHz, with the Caltech SubmiUimeter Observatory (CSO) in 
January 2011. The beam size was 32" FWHM and the beam efficiency was 0.73. The full 
spectrum is shown in Figure 4. We have focused on the relatively unblended and unconfused 
trio of hyperfine components near 250.8 GHz, and used CLASS software to fit the spectrum 
based on optically thin relative line intensities and two line components. The parameters 
for the wide component arc vq = 9.3 km s^"*^, AfpwHM = 19.0 km s~-^, and for the narrow 
component, Vq = 10.4 km s~^, and AwpwHM = 3.1 km s~^ The data and the best-fit profile 
with these two components are shown in Figure 6. The excellent fit confirms that there is a 
component having centroid velocity and line width quite different from those characteristic 
of the Hot Core or the region around KL, and which are close to the parameters of the O2 
fines. 

The integrated intensity in the narrow component of the NO 250.796 GHz line is 
jTyiBdv = 1.4 K km s~^ Non-LTE models of NO were carried out using the RADEX 
code and collision rates from Lique et al. (2009b). They indicate beam-averaged column 
densities of 2.7 x 10^^ and 1.45 x 10^^ cm-^ for ^(Ha) = 5 x lO'^ cm'^, Tk = 77.5 K, and 10^ 
cm^^, 100 K, respectively. The corresponding column densities of O2 in a non-LTE treat- 
ment are 6.5 and 7.5 xlO^® cm~^, respectively. Apparent beam-averaged abundance ratios 
N0/02= 0.04 and 0.02 arc deduced. The excitation of NO is more sensitive to density and 
temperature than that of O2. Observations of higher transitions of NO will help to constrain 
the conditions in the 02-emitting region and ultimately could help determine the ratio of 
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Fig. 5. — Spectrum of O J = 7 - 6 emission at the H2 Peak 1 position obtained simul- 
taneously with the O2 774 GHz data using Herschel HlFl. The four-Gaussian fit indicates 
a pedestal feature ~ 22 km s~^ wide (long-dashed curve) and three narrower components. 
The narrowest of these components peaks at a velocity of 10.3 km s~^. 
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atomic nitrogen and oxygen in the gas phase. 

HIFI observations of the 620.710 GHz 832 transition of ortho-H20 in Orion have 
yielded, in addition to a broad thermal emission feature seen throughout the central portion 
of the cloud, a very narrow feature close to 12 km s~^(Melnick et al. 2011, in preparation). 
This has many characteristics of maser emission, with extremely narrow lineswidths < 1 km 
s~^. The size of the narrow-line emission region is difficult to determine due to the presence 
of large gradients in the thermal emission, but the source location is likely coincident with or 
slightly north of the H2 Peak 1 position. Populating the upper level of this transition, more 
than 730 K above the ground state (as well as the location of the emission) is more consistent 
with spatially-extended shocked gas than with any of the compact sources discussed above. 

3.3. Constraints on Source Location From Pointing Offsets 

The two linear polarizations (H and V) observed by HlFl have a small relative pointing 
offset, which is different for each band. There are clear detections in each polarization of 
the 774 GHz O2 transition each of two well-separated observing epochs (March and Septem- 
ber 2010). The 1121 GHz line lacks sufficient signal to noise ratio to allow epoch- and 
polarization-dependent analysis, and the 487 GHz transition has a much larger beam and 
marginal signal to noise ratio. 

When processed through HIPE version 5, the pointing offsets of the polarizations rel- 
ative to the nominal pointing direction are recorded. For March, the a and S offsets for 
H and V are (-0.9",+2.1") and (+0.9",-2.1"), while for September they are (+0.5",-2.2") 
and (-0.5", -1-2. 2"), respectively. Thus, the offsets from separate polarizations have almost 
been exactly interchanged due to the 6 months that have elapsed between the two sets of 
observations. 

The V(March) and H(Sept) pointing direction is in a direction offset towards the Hot 
Core by 2.3", and the V(Sept) and H(March) data are offset in the opposite direction by 
an equal amount. We can investigate whether the measured intensities are consistent with 
a source for the O2 emission located at the position of the Hot Core. To do this, it is thus 
reasonable to average together the H data from September with the V data from March, and 
the V data from September with the H data from March. Doing this we find {J T[H(Sept) 
+ V{Mar)]dv) = 0.160 ± 0.023 K kms^^and (/ T[V(Sept) + H(Mar)]dv) = 0.210 ± 0.012 K 
kms-^ The ratio of the integrated intensities is (/ T[V(Sept), H(Mar)]dv)/(/ T[H(Sept), V(Mar)](iv) 
= 1.31± 0.20. 

The Hot Core is separated from the H2 Peak 1 nominal pointing position by 24.6". 
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Fig. 6.— Spectrum of the F = 7/2 - 5/2, F = 5/2 - 3/2, and F = 3/2 - 3/2 hyperfine 
components of the J = 5/2 - 3/2 transition of NO. The relative intensities and offsets 
referenced to the best-fit velocity of 10.4 km s~^ for the low dispersion portion of the strongest 
F = 7/2-5/2 hyperfine component are indicated by the vertical lines. The combined best-fit 
spectrum of the two kinematic components for each of the three hyperfine components is 
indicated by the solid green curve. 
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Adding the offsets for the different epochs appropriately, we find that the distance for the 
H(Sept) and V(Mar) data is 22.4", and for the V(Sept) and H(Mar) data is 26.9". Taking 
the Hot Core Source to be small compared to the antenna beam size, the ratio of the signal 
measured when pointing closer to the source to that found when pointing farther from the 
source is given by exp{{dfarther/OoY — {dcioser/OoY), where 9o is the 1/e radius of the beam, 
equal to 16.6" at 774 GHz. The resulting ratio is 2.27. 

Relevant information about pointing and source location is provided by the other lines 
observed simultaneously with the O2. We see changes in most lines, and the C^^O line is 
particularly useful as it is relatively isolated as well as quite strong. It seems reasonable that 
the broad component of this line originates in a relatively small region near the Hot Core. 
Table 3 includes the data for the two polarizations at two epochs. The strength of the broad 
component is sufficient that we find ratios of the integrated intensities (as defined above) 
V(Mar)/H(Mar) = 1.41 and H(Sept)/V(Sept) = 1.77, with an average value equal to 1.59. 
Both ratios are in the sense expected for a point source in the vicinity of the Hot Core, and 
"flip" in the manner that is expected. 

The formal statistical error for each ratio above is less than 4 percent, so these two 
measurements are marginally consistent with their average value. However, the data in 

Table 3 suggest that there is another systematic difference, which can be seen in the average 
of H and V in September being significantly greater than the average of H and V in March 
by a factor of 1.5. Since the relative pointing offsets of the two polarizations arc essentially 
just interchanged between the two observing epochs, this suggests that there was an overall 
pointing change between the periods in which the observations were carried out. 

Assuming we are observing a pointlike source at the nominal distance of the Hot Core, 
we can solve for the pointing offset, which is 2.4" in the sense that the average pointing 
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direction in September was this amount closer to the position of the Hot Core than in 
March. Although this error results in a large difference in the signal received when observing 
a source on the "side" of the beam, this is a relatively small pointing error (~ FWHM/12), 
which is within the expected pointing precision for Herschel (Pilbratt et al. 2010). Our data 
are only marginally consistent with the O2 emission being a small source at the location of 
the Hot Core. The NH3 emission from the northwest portion of the Hot Core, and the Peak 
A - Western Clump, and the Murata et al. CS clumps are all somewhat closer, and thus 
more consistent with our measurements. 



4. Analysis 

4.1. Collisional Excitation of O2 

Since the millimeter- and submillimeter- wavelength transitions of O2 are magnetic dipole 
transitions, their spontaneous decay rates are relatively small. Thus, in the relatively dense 
regions being observed, it is likely that collisions are the dominant excitation mechanism. 
Orlikowski (1985) carried out calculations for collisions between O2 and He atoms (as sur- 
rogates for H2 in J = level) for levels defined by rotational quantum numbers N = 1 and 
3, at a total energy of 300 cm~^ (corresponding to 432 K). Corey, Alexander, & Schacfer 
(1986) performed similar calculations for excitation from N = 1 and 3 to N = 3, 5, and 
7 levels for a collision energy of 313 K. The results suggest a propensity for collisions in 
which the quantum number F, which is equal to N - J (and hence = -1, 0, or -|-1) does not 
change. Bergman (1995) and Valkenberg (1995) carried out close-coupled calculations of 
O2 by He over a range of temperatures relevant for interstellar clouds, using a semi-empirical 
interaction potential and a number of approximations. 

The uncertainty due to the potential energy surface (PES) has been reduced by use of 
an improved surface in the recent close-coupled calculations by Lique (2010). The recent 
study finds that the collisions in which F is unchanged have rate coefficients ~10~^° cm^ s~^, 
essentially independent of N. The collisions in which N changes are a factor of 10 smaller 
for low N, and almost a factor of 100 smaller for N > 10. The transitions with the largest 
cross sections do not correspond to ones for which radiative transitions are allowed. Thus, 
the upper level of an observed transition is generally not populated by a collision from the 
lower level of the transition. The result is that a critical density defined in terms of the 
spontaneous decay rate divided by the deexcitation rate coefficient (for the same levels) is 
overestimated. 

Given that we expect the transitions observed to be optically thin, the most meaningful 
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quantity is simply the fractional population of the upper level, since the observed antenna 
temperature is directly proportional to this quantity. Using the rate coefficients at a kinetic 
temperature of 100 K from Lique (2010), scaled by a factor of 1.37 to account approximately 
for H2 rather than He as a collision partner, we find that the upper level fractional populations 
are within 10% of their LTE values for densities above Ix 10^ cm~^ for the 487 GHz transition, 
1.3x10^ cm~^ for the 774 GHz transition, and 2x10"^ cm~^ for the 1121 GHz transition. This 
density for the 1121 GHz transition is a factor of 6 smaller than the critical density defined 
in the usual manner. 

In moderate density interstellar clouds, the molecular hydrogen density is n(H2) — 10^ 
cm~^, which is sufficient to bring only the population of the (3,3) level close to LTE; the 
populations of the upper levels of higher-frequency transitions will be subthermal. However, 
in regions such as the centers of giant molecular clouds, where densities are generally > 10^ 
cm~^ and certainly in excess of 10^ cm~^, the fractional populations of the upper levels of 
all transitions will be close to LTE values. This is confirmed by a fully non-LTE treatment 
of the O2 excitation, which produces the same result for column density and temperature 
when the hydrogen density is > 10^ cm^^. It is thus appropriate to use the relative line 
intensities with the assumption of LTE to determine the kinetic temperature and the total 
O2 column density. 



We can use the three transitions observed to determine the temperature of the O2- 
emitting region. To do this we first assume that the source can be characterized by a single 
temperature, and that it fills the beam of all three transitions observed. Both of these 
assumptions may only hold approximately. If there is a fractional beam filling factor less 
than unity, but which is the same for all beams, it does not affect their relative intensities 
or the beam-averaged column density that we determine. Assuming that the source does fill 
the main beam of all three transitions and that the emission is optically thin, the integrated 
main beam brightness temperature is proportional to the upper level column density of the 
transition observed, A^„, through the relationship 



where v is the frequency of the transition and A^i is its spontaneous decay rate. Assuming 
LTE, as discussed above, the ratio of the column densities of two transitions is determined 
only by the kinetic temperature of the gas. With three transitions, we can form two such 
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ratios, and from each determine the range of kinetic temperatures consistent with the ob- 
servations (including at least statistical uncertainties in J Tmbdv). The results arc shown in 
Figure 7. The data in Tabic 1 yield 0.45 < I(487)/I(774) < 0.63 and 0.40 < I(1121)/I(774) 
< 0.65. These limits define the temperature limits shown in Figure 7. A kinetic temperature 
between 80 K and 100 K seems most likely, but values between 65 K and 120 K are consistent 
with the relative intensities of all three transitions observed. 

From equation 1, we obtain that N(5,4) = 4.5 ± 0.21 x 10^^ cm~^. With the kinetic 
temperature defined consistently by the two ratios of the three integrated intensities, the 
total O2 column density relative to that in (N, J) = (5,4) can be reliably calculated, and 
the result is shown in Figure 7. The energy of the (N, J) = 5, 4 level is such that there is 
remarkably little variation in the fraction of the total O2 in this one level. We combine the 
uncertainty in N(5,4) with that in the kinetic temperature to determine the beam-averaged 
column density N(02) = 6.5 ± 1.0 xlO^^ cm^^. 

4.3. Filled Beam Fractional Abundance of O2 

Due to the wide range of temperatures and kinematic variations in the gas towards H2 
Peak 1, it is difficult to determine the H2 column density to compare with that of the O2 
observed. The C^^O J = 7 - 6 transition discussed in §3.2 is a good candidate for this, since 
it is optically thin, and at these elevated temperatures, depletion of carbon monoxide should 
not be significant. The width of this line is sufficient to include all hyperfine components, 
but as discussed in §3.2, we do see three velocity components in addition to the very wide 
component, presumably associated with the molecular outflow. A four- Gaussian fit is very 
satisfactory, and yields the results given in Table 2. While the velocity of component 4 (10.3 
km s~^) agrees best with that of the O2, the line width (1.3 km s~^) is considerably narrower 
than the ~ 3 km s~^ measured for the three O2 lines (Table 1), so it is difficult to make a 
definitive correspondence. 

Using equation 1, wc find J = 7 C^''0 column densities 3.3x10^^ cm~^ for the 10.3 km 
s~^ component and 7.8 x 10^^ cm^^ for the combination of the three narrow components (2 to 
4). If we consider the emission at velocities above 10 km s~^ excluding the broad components 
1 and 2, we obtain a column density equal to 2.5 xlO^^ cm~^. The fraction of C^^O in J 
= 7 becomes almost independent of temperature for T > 80 K (varying from 1/13 at 80 K 
to a minimum of 1/10 at 150 K, and increasing very slightly higher at temperatures below 
300 K). Wc adopt a value of 1/11, which yields N(C^^O) summed over all levels equal to 
3.6x 10^^ cm^^ for the very narrow component and 8.7x 10^^ cm~^ for the narrow components 
combined. We take the ratio CO/C^^O = 2460 (Penzias 1981), and H2/CO = 10^ yielding 
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Fig. 7. — Lower panel: Intensity of the 487 GHz O2 transition relative to the 774 GHz 
transition, and of the 1121 GHz transition relative to the 774 GHz transition. These values 
assume that the beam filling factor is the same for the three transitions, and that all levels 
are populated according to LTE. The triangles are from the results given in Table 1. The 
dotted boxes are defined by the statistical uncertainties in the observed intensity ratios, 
and the corresponding values of kinetic temperature that are implied. Upper panel: Total 
column density of O2 relative to that in the upper level of the 774 GHz transition, (N,J) = 
(5,4). In LTE, this depends only on the kinetic temperature as shown here. 
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N(H2) = 8.9x1021 cm-2, G.TxlO^^ cm-^ and 2.1x1023 crn-^, for the 10.3 km feature, the 
>10 km emission, and for the combined narrow components, respectively. 

This range of column densities corresponds to a range of densities averaged over the 
projected diameter of the largest Herschel beam (44") of (n(H2)) = N{E.2)/GD = 3 x 
10^ to 8 X 10^ cm~3. The distance D — 416 pc is the weighted mean of four recent parallax 
measurements of masers in the Orion KL region and of stars in the Orion Nebula Cluster 
(Hirota et al. 2007; Sandstrom et al. 2007; Menten et al. 2007; Kim et al. 2008). These 
densities are in harmony with the requirements of the preceding LTE analysis of the O2 
line intensities. With these values bracketing the range of H2 column density associated 
with the O2 emission, we find that the fractional abundance of O2 is between 3.1x10"^ and 
7.3x10"^. The higher value is still only ~ 1% of the total oxygen available. The lower value 
is comparable to the upper limits found by SWAS and Odin, if the fractional abundance of 
carbon monoxide is lower (which is plausible, the adopted value being close to the maximum 
observed), the column density of H2 will be increased and the fractional abundance of O2 
proportionately reduced and thus be even smaller compared to the total elemental oxygen 
abundance. 



5. Beam Filling and Source Offset Effects 

Modeling of the emission for anything but an extended source with a single filling factor 
for all transitions involves source coupling to beams of the three transitions each having a 
different beam size. This leads to some effects which we discuss here. 



5.1. Effect of Beam Filling 

A beam filling factor less than unity, but which is the same for all transitions, can be 
produced by an extended but clumpy medium (as suggested by e.g. the shocked H2 emission 
seen in Figure 2). This situation leaves the relative intensities of the transitions unchanged, 
and thus the temperature we determine for the region responsible for the O2 emission will 
be unaffected. The column density of O2 in the "filled" portion of the beam is increased by 
the reciprocal of the beam filling factor. If the O2 emission is associated with postshock gas 
traced by H2 emission, then as shown in Figure 2, the emitting region is extended on the 
scale of all of the HIFI beams, albeit modestly weaker in the outer portions of the lowest 
frequency (487 GHz) beam. The fractional abundance of O2 may be increased if the tracer 
of H2 is uniformly distributed, but may be essentially that derived above if the overall mass 
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distribution is nonuniform, rather than only the distribution of O2. 

If Peak A is the source of the observed O2 emission, we have to deal with a small source 
which is also away from the antenna boresight direction. In Appendix 2 we discuss the 
limiting case of a source that is so small that it can be considered to be point-like relative 
to all beams. The ratio of the intensities of two transitions, 2 and 1, can be considered to 
be the ratio characteristic of the source itself multiplied by the correction factor R(2/l) as 
defined in equation A2-7 in Appendix 2. To derive the required characteristics of the source, 
we must thus divide the observed ratio by the correction factor. For the ratio of the 1121 
GHz transition to the 774 GHz transition we find 

/°^-(1121) _ /-(1121) . A(774) _ ^ ^ /-(1121) 
/ofe(^774) prc^77A) ^A(1121)^ ' P^'^(774) ' ^ ' 

The ratio for the source has to be a factor ~ 2 smaller than the observed ratio, which as seen 
in Figure 7, means a lower kinetic temperature. For an observed ratio of 0.5, a temperature 
of 70 K is required for a source filling both beams. However, the ratio for a small source 
must be smaller by a factor of 2.1, and hence equal to 0.24. As seen from Figure 7, this 
requires a kinetic temperature of ~ 35 K. 

The second hue ratio is again normahzed to the intensity of the 774 GHz transition, 
and hence the correction factor it:(487/774) = [A(774)/A(487)]^ = 0.4. As indicated by the 
appropriate curve in Figure 7, this again corresponds to a lower kinetic temperature for the 
small source. An observed ratio of 0.5, which requires a filled source temperature of 100 

K demands a ratio 2.5 times larger in a small source, and thus a source temperature again 
~ 30 K. A small source on boresight is consistent with a single temperature determined by 
the ratios of integrated intensities, but the source temperature is too low given the case of 
excitation of O2 and the temperature of the center of the Orion molecular cloud. It is thus 
unlikely that the observed O2 lines are produced by a small source coincident with Orion H2 
Peak 1, but an extended while beam-diluted source can fit the present data. 



5.2. Effect of Source Offset 

We again restrict our discussion to a small (point-like) source, which when offset by an 
angle 9 produces a signal reduced from that of a source on boresight by a factor which is 
given by the normalized antenna response pattern 



Pn = e-^'/'°^' , (3) 
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where (= 0.6019fwhm) is the 1/e beam radius. 9o is proportional to the observing 
wavelength, and so a given offset produces a greater reduction in observed integrated intensity 
at a higher frequency. Following the discussion of beam filling, we define the integrated main 
beam temperatures with the source off and on boresight as 1°^ and /°", respectively. We 
then find 

7°g(1121) _ eWMm? 7°n(ii21) 

Since ^^o(1121) = 11.4" and 6o{774) = 16.6" the correction factor for off-borcsight point- 
ing is always less than unity and decreases as the offset increases. To fit a given observed 
ratio, we require an on-boresight source with a larger ratio, which can only be produced by a 
higher source temperature. The ratio for the 487 GHz transition to the 774 GHz transition 
shows the opposite dependence, being greater than unity. Thus, to fit an observed integrated 
intensity ratio, the off-boresight source would need to have a smaller 487 GHz to 774 GHz 
ratio than that for an on-axis source, and the effect again is that the source temperature 
would be greater than if the source were on boresight. 

The observed 1121 to 774 ratio, close to 0.5, is what would be produced by a on-boresight 
small source that if filling both beams would produce a ratio equal to 0.5/2.1 = 0.24. To 

obtain a ratio for a small source off-boresight equal to the on-boresight ratio for a source that 
fills the beams requires that the off-borseight correction be the inverse of the small-source 
correction, and thus equal to 1/2.1 = 0.48, which occurs if the source is off axis by just over 
13". 

As seen in Figure 7, the ratio J^^'^(1121)/J*^'^(774) has an upper limit of about 0.8 for a 
filled source, irrespective of temperature. This would be produced by a on-boresight small 
source with an (unphysical) ratio of 1.68. However, a very warm small source off'-boresight 
could produce the observed ratio of 0.5 with off-boresight correction factor equal to 0.30, 
which results when the source is off boresight by about 17". Especially when the uncertainties 
are considered, consistent solutions for all three transitions are possible for small sources over 
a restricted range of offset angles, up to approximately this hmit. Since the offset to the 
Hot Core is close to 27 it appears difficult to have this source dominate the O2 emission 
we have observed at H2 Peak 1. One or more of the compact sources discussed previously 
are more likely the source of the emission. In particular. Peak A (with an offset of 21"), is 
consistent with this picture, with gas and dust temperatures well over 100 K being indicated 
by other observational data and required to fit the present observations of O2. 
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6. Discussion 

There seems little reason to doubt that we have detected emission from O2 but the source 
of the emission remains imperfectly characterized. The issue at hand is to understand the 
circumstances that result in the increase in the O2 abundance from the low levels derived from 
SWAS and Odin observations. Two theories for enhanced abundance of O2, (1) warm dust 
with restoration of pure gas-phase chemistry, and (2) shocks, seem to be relevant possibilities 
for explaining the present observations. We here give some additional details that bear on 
explaining our data for H2 Peak 1 in Orion. 

6.1. Warm Dust 

It is generally accepted that dust grains in dense clouds without star formation in 
the immediate vicinity will be sufficiently cold that over time (depending on the density), 
atomic oxygen will deplete onto grains and to a large degree be hydrogenated to water. In 
consequence, the grains acquire water ice mantles incorporating a significant fraction of the 
available oxygen, with an accompanying significant reduction in the gas-phase X(02), as 
developed in the models by Bergin et al. (2000) and Hollenbach et al. (2009) discussed in 
§1. The formation of a massive star in or near the region in question produces radiation that 
is in a short distance degraded to infrared wavelengths, and which heats dust grains. At 
modest temperatures, any desorbed O2 will be returned to the gas, as discused in §1.1. When 
the grain temperature reaches ~ 100 K, return of water ice to the gas phase becomes rapid 
(Eraser et al. 2001). This starts the restoration of "pure" gas-phase oxygen chemistry. The 
time required to reach steady state depends primarily on the ionization rate, but ultimately 
a very substantial X(02) > 10~^ will be established. 

In Figure 8, we show results of a representative calculation, which is based on NAHOON, 
the Wakelam et al. (2005) model, but using the OSU-1-2009 reaction rate file. The H2 
density is 10^ cm~^, the gas temperature is 100 K, and grain surface reactions are excluded 
except for H2 formation. We assume that at t = 0, all of the elemental oxygen that had 
been on the grain surface is returned to the gas phase, and that the oxygen that had (in the 
steady-state solution for large times in this model) been in molecular and atomic form, as 
well as that having been water ice, comes off the grain as H2O. This reflects the idea that 
at some point in the past, the grains had been cool enough that all of these species would 
deplete onto the grain surfaces, and end up as water ice, but the result here is in effect a 
gas-phase chemistry model. There are also trace oxygen-bearing species such as SO and SO2 
that are in the gas phase at the start of the present calculation, but together they represent 
less than 0.5% of the total available oxygen. 
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Starting at t = 0, the O2 fractional abundance builds up steadily, and although initially 
smaller than that of atomic oxygen, soon surpasses it. The two species end up with com- 
parable steady-state abundances 5 to 7 xlO~^. The time to achieve maximum abundance 
of O2 scales inversely with the cosmic ray ionization rate. The time scales for achieving a 
given O2 fractional abundance as well as steady-state abundance are largely independent of 
the density of the region. The behavior is also largely independent of temperature as long 
as the dust is warm enough to prevent any significant grain surface depletion to occur. The 
steady-state fractional abundance of H2O increases as the ionization rate increases, although 
the abundances of O and O2 are only slightly affected. 



Fig. 8. — Time evolution of gas-phase abundances of atomic oxygen, molecular oxygen, 
water, and carbon monoxide in a cloud with dust and gas temperatures equal to 100 K and 
H2 density equal to 10^ cm~^. At t = 0, elemental oxygen not in CO is almost entirely in 
the form of gas-phase water, grain surface water having been returned to the gas phase as 
a result of the warm grain temperature. The dotted curves are for a cosmic ray ionization 
rate of 10~^^ s~^, the dashed curves for 10^^^ s^^, and the solid curves for 10~^^ s~^ 

To evaluate the quantity of gas that may be maintained at conditions for the above 
chemistry to evolve, we evaluate a highly idealized model based on a single luminous source 
in a spherically symmetric cloud. Taking a central source with 10^ - 10^ Lq, plausibly IRc2 
(Downes et al. 1981) or else a collection of sources in the region (Gezari, Backman, & 
Werner 1998) having similar total luminosity, we have carried out a representative cal- 
culation of the dust temperature distribution using the DUSTY radiative transfer code 
(http://www.pa.uky.edu/~moshe/dusty; Nenkova et al. 2000). The center to edge visual 
extinction is equal to 1000 mag., but this value is unimportant as long as the dust at the 
cloud edge is cooler than the temperatures of interest. We find that the solutions at distances 
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of interest are quite independent of the radial density distribution assumed to be a power 

law n{r) ~ r~" with index a between and 2. At the 4x10^^ cm projected distance of 
Peak A from lRc2, the dust temperature is 56 K for the lower luminosity and 100 K for the 
higher. At the 1.5x10^'' cm projected distance of H2 Peak 1 from lRc2, the temperatures 
are 35 K and 55 K, respectively Thus, even if Peak A is externally heated, the grains will 
be sufficiently heated to result in "warm dust" chemistry, while those at the distance of H2 
Peak 1 will not. 

The fraction of the total opacity between the central source and the location of a partic- 
ular grain temperature docs depend on a. For relatively fiat density distributions (a < 1), 
the fraction of the total optical depth is < 0.5 for a dust grain temperature of 100 K. While 
recognizing that this idealized, spherically symmetric model may not be accurate for this 
complex region, it seems plausible that IRc2 or another source, or the combination of sources 
in the region can result in warm grains in Peak A, even if this source is not self-luminous. 
The greater distance to H2 Peak 1 means that the grains there would still be cold enough 
to retain significant ice mantles if IR heating were the only operative process warming the 
dust. 

Given that the grains in Peak A may be sufficiently warm that the gas-phase O2 abun- 
dance has been significantly raised, we can examine whether this source is consistent with 
the O2 line intensities we have observed. We must thus consider a small source offset from 
the beam pointing direction. Following the discussion in §5.1 and §5.2, the 1121/774 ratio 
is increased by a factor of 2.1 due to small source size and decreased by factor 0.15 due 
to offset, relative to fiUed-beam on-boresight model. The resulting correction factor is 0.3, 
which means that the maximum ratio for a warm beam-filling source of 0.8 becomes 0.25, 
lower than the observed ratio of 0.51±0.1. For the 487/774 ratio, we get a combined cor- 
rection factor of 1.12, which when multiplied by the high temperature filled and centered 
model ratio of 0.45, gives a ratio equal to 0.51, compared to the observed ratio of 0.54±0.1. 
Considering the uncertainties in source size and geometry, in addition to calibration issues, 
we conclude that the observed line ratios are consistent with a small warm source located at 
the position of Peak A. 

If Peak A is the source of the 11 km s~^ to 12 km s~^ emission, the column densities 
are certainly greater than the filled beam values given previously. With a representative 10" 
for the diameter of Peak A (larger than Murata et al. (1992), as there is some extended 
emission) and including the attenuation due to being 21" away from boresight, we obtain a 
total correction factor of 55 compared to a fiUed-beam on-boresight source. The C^^O used 
to determine the associated H2 column density certainly does not all come from this small 
source, but scaling the H2 column density impfied by the > 10 km s~^ emission, 6.7x10^^ 
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cm~^, by the correction factor, we find N(H2) = 3.7x10^^ cm~^. The conversion from C^^O 
in J = 7 to total N(C^^O) is essentially independent of temperature for T > 100 K. The 
origin of the O2 emission in this small region leaves the emission optically thin, as is the 
C^^O emission for v > 10 km s~^. For a source temperature > 150 K compared to the ~ 
100 K value in the filled-beam model, the column density of O2 will be slightly increased 
due to the larger partition function (Figure 7). For a temperature of 180 K we obtain N(02) 
= 4.6x10^^ cm~^ for a 10" diameter source. The fractional abundance of O2 is only shghtly 
changed and is equal to 1.3x10"^, a value consistent with the warm grain chemical models. 
According to the warm dust model as shown in Figure 8, for a low ionization rate of 10~^^ 
s~^, the time to reach this fractional abundance is < 10^ yr. 

Murata et al. (1992) derive a mass of 3.2 M0 for Cut D clump, coincident with the 
Peak A - Western Clump - MF4 source, as discussed above. Wright et al. (1992) carried 
out submillimeter dust observations with somewhat lower angular resolution, and derived a 
size of 5" ~ 10 " and a mass of 4 - 13 M©, slightly larger than those of Murata et al. (1992), 
possibly due to their lower angular resolution. The peak column density from the Murata et 
al. (1992) data is N(H2) = 5x10^^ cm~^, which agrees surprisingly well with that implied by 
our C^^O and related observations keeping in mind that the mass derived from submillimeter 
dust emission has its own ~ factor of three uncertainty. The average density in the clump 
is ~ 2x10^ cm~^, which is certainly adequate to thermalize both the O2 and the C^^O. 

The simplest estimate for the virial mass of this clump, assuming uniform density, and 
ignoring magnetic fields and external pressure, is M^i- = 5af^^K/G, where aios is the line of 
sight velocity dispersion and R is the cloud radius. The result is M^j^ = 8-19 Mq, and 
Myir/M = 0.7 - 4 for the range of sizes and masses determined from the dust emission. This 
range includes values characteristic of cores in general. While we cannot definitively prove 
that the observed O2 lines arise from this compact source, it certainly does seem to be a 
good candidate to explain our observations. 

6.2. Shocks 

The H2 Peak 1 position was selected as it coincides with the strong local maximum in 
the intensity of the v = 1 - S(l) transition of molecular hydrogen studied by Beckwith 
et al. (1978). Using several of the vibration-rotation transitions, these authors found an 
excitation temperature in excess of 2000 K. Since the peak is not in the immediate vicinity 
of any luminous source, this very high value is a strong indication that this is a result of 
shock excitation. In the intervening years there has been a plethora of papers that report 
dramatically improved spatial and spectral resolution (Scoville et al. 1982; Chrysostomou 
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et al. 1997; Schild, Miller, & Tennyson 1997; Stolovy et al. 1998; Schultz et al. 1999; Salas 
et al. 1999; Rosenthal, Bertoldi, & Drapatz 2000; Vannier et al. 2001; Gustafsson et al. 
2003; Kristensen et al. 2003; Lacombe et al. 2004; Bally et al. 2011). 

Although fluorescence powered by UV radiation from the Trapezium stars dominates 
the H2 line emission on degree scales in Orion (Luhman et al. 1994), the emission of highest 
surface brightness in the molecular core is due mainly to coUisional excitation in shocks 
(Rosenthal, Bertoldi, & Drapatz 2000). As discussed by these authors, population diagrams 
based on the wide range of lines available from ISO clearly indicate a large range of excitation 
temperatures spanning a few hundred to several thousand K. Multiple shocks are indicated 
by the pure rotation H2 transitions observed by Parmar, Lacy, & Achtermann (1994). A 
relatively low-temperature component is also indicated by the 17 /im J = 3 - 1 rotational 
transition studied by Burton & Haas (1997). This study also suggested the presence of 
shocks having velocities as low as 5 km s~^. 

High resolution observations provide evidence of a complex kinematic structure. Among 
the highest spectral resolution data are those of Chrysostomou et al. (1997) of the H2 v = 1 
- S(l) line, carried out with a velocity resolution of 14 km s~^. These data indicate that the 
spatially extended H2 emission peaks at a velocity of 12 km s~^ and has an average FWHM 
line width equal to 37 km s~^, but that localized regions have different central velocities and 
broader line widths. The highest spatial resolution observations (0.15", Gustafsson et al. 
2003) show complex localized structures having dramatically different central velocities. 

The effects of shocks on the abundance of O2 have been considered in some detail by 
Bergin, Melnick, & Neufeld (1998); references to reaction rates and earlier studies can be 
found in that paper. The key point is that the reaction O + H2 — OH + H, which has 
an activation barrier of 3160 K, can become the dominant production pathway for OH at 
postshock gas temperatures. This provides a major additional source of OH beyond the 
primary pathway operative in quiescent clouds in which dissociative recombination of H3O"'" 
is the major source of OH. Thus, if there is appreciable preshock atomic oxygen, there will be 
enhanced production of OH in the shock, which then can react with O via the temperature- 
independent reaction OH -|- O — > O2 -l- H discussed in §1. FIR observations (Goicoechea 
et al. 2006) suggest that there is large increase in the abundance of OH as a result of the 
shock in the Orion KL outflow, but the velocity and line width of the OH do not allow any 
conclusive association with the O2 emission. 

However, the reaction OH -|- H2 — )■ H2O -|- H, with an activation barrier of 1660 K, 

competes for the available OH. In addition, there is a back reaction H + O2 — OH + O, 
which limits the abundance of O2 if the gas temperature exceeds ~ lOOOK and there is 
some atomic hydrogen present. These conditions are likely to occur following the passage 
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of a strong nondissociative shock. Thus, shocks strong enough to enable production of OH, 
but in which the postshock temperature does not exceed ~ 500 K, result in maximum O2 
production (Bergin, Melnick, & Neufeld 1998). Kaufman (2010) has studied a wide range 
of nondissociative C-shocks over a wide variety of shock velocities and preshock densities, 
and finds that maximum O2 production results from shock velocities 10 km s^^ to 12 km 
s~^, and is relatively independent of preshock gas density. The range of shock velocities over 
which O2 column densities in excess of 10^^ cm~^ are produced is broadened for n(H2) in the 
preshock gas > 10^ cm~^. Under optimal conditions, O2 column densities can be as large as 
10^^ cm-2. 

The models of shock chemistry, while still highly idealized relative to the complex struc- 
ture in the H2 Peak 1 region (see e.g. Figure 2), may well produce a relatively large column 
density of O2. The conditions required are (i) a C-shock of appropriate velocity, (ii) low 
electron density with X(e~) < 10~^ - 10~^, (iii) sufficiently high preshock density, n(H2) > 
10^ cm~^, and (iv) a significant fraction of available oxygen in gas-phase atomic form. The 
first three requirements do not seem implausible in this region of Orion, due to the high de- 
gree of dumpiness, the high average density, and the wide range of shock velocities present. 
In considering shock models, we have to allow for a range of shock velocities propagating 
in a very nonuniform medium. This is to some degree the picture advanced by Salas et al. 
(1999) in which shocks driven by a stellar wind impinge on (possibly moving) clumps. The 
shock velocity in a dense clump will be significantly lower than that in the diffuse interclump 
medium, and so one can easily imagine a wide range of shock velocities. In addition to the 
shock propagating into the clump (which may have too low a velocity to produce abundant 
O2), there can be a variety of other complex shocks that could release atoms and molecules 
from the grains, which could subsequently form large amounts of O2. 

The issue of the atomic oxygen abundance is perhaps the most challenging one. For 
C-shock speeds in excess of ~ 15 km s~^, the ion-neutral streaming is sufficient to sputter 
grain mantles and thus produce a large density of atomic oxygen. Even in the absence 
of sputtering, if the extinction of this region relative for the UV field produced by the 
Trapezium stars is not greater than ~ 9 mag, the atomic oxygen fractional abundance could 
exceed 10^^, which would be sufficient for the C-shock to produce observed column density 
of O2. The O2 column density in the PDR itseff (< 10^^ cm-=^; HoUenbach et al. 2009) 
is significantly below that observed unless there is some very special geometry involved. 
Alternatively, a higher-velocity J-shock (such as suggested by the H2 data) could produce 
the atomic oxygen. In the dense postshock gas dissociated H2 reforms quickly. A lower- 
velocity C-shock propagating into this gas would then produce the O2 as discussed above. 
It is also possible that UV radiation generated in the J-shock can, through photo-reactions 
induced in water ice molecules on grains, result in ~ 10% of the molecules desorbed being O2 
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(Oberg, van Dishoeck, & Linnartz (2009) and Oberg, private communication 2011), which 
could be a supplement to the other production channels. 

The hne profile of the O2 produced in a shock depends in detail on the geometry, but 
for the ideahzed case of a planar low velocity shock, has a quite narrow hne "core" together 
with a broad "wing" profile. In order for the hne velocity to be shifted from that of the 
ambient material in the center of Orion by only 3 to 5 km s""*^, the shock direction would 
have to be close to the plane of the sky. However, if an initially planar shock encounters 
inhomogeneities, the propagation direction will be modified, and it seems inevitable that 
some line broadening would result. The observed 11 - 12 km s~^ velocity is shifted relative 
to the 8-9 km velocity of the ambient gas in the opposite sense of the CO outflow, 
which has its blue-shifted wing extending towards H2 Peak 1 (Erickson et al. 1982). 

A serious problem with the shock model is that while the observed N(02), 6x10^^ 
cm~^, is ~ a factor of two less than the maximum N(02) that can be produced, the observed 
column density is a beam-averaged value. Thus, we require the shock to essentially fill a ~ 
30" region, corresponding to ~ 2x10^^ cm in size at the 416 pc distance of Orion. Given 
the highly inhomogeneous structure of the H2 emission seen in Figure 2, it is difficult to 
see how our three HIFI beams could be filled with gas that has been through a low- velocity 
shock all giving projected velocities within a ±1.5 km s~^ range. It is possible that the 
observed emission could be the result of shocks having passed through beam-diluted clumps 
distributed throughout the beam, thus allowing a reasonable agreement with the relative 
intensities we have observed. The O2 column density would, however, have to be higher by 
the inverse of the beam-filling factor. If this factor were as large as few, it would put N(02) 
close to the upper limit for shock production of this species. 

In summary, while shock models can in principle produce the observed column density 
of O2, there are difficulties regarding the initial conditions, particularly the abundance of 
atomic oxygen, and the observed narrow line profiles. These can be addressed only by more 
detailed and extensive modeling. 

7. Conclusions 

We have observed three submillimeter rotational transitions of the O2 molecule in emis- 
sion, using the HIFI instrument on the Herschel Space Observatory. The pointing position 
was towards Peak 1 of the H2 emission in Orion. Signals were detected from the 487 GHz, 
774 GHz, and 1121 GHz transitions, having peak velocities between 11 km s~^ and 12 km 
s~^, and FWHM hne widths ~ 3 km s~^. The detection of three transitions and the ehmi- 
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nation of the most obvious possible "interloping" lines confirm that the emission is from O2. 
The determination of the source of the emission and its characteristics are difficult due to 
the complexity of the region, but the highlights of the analysis are the following. 

1. Assuming that the emission fills the beam at all three frequencies, the relative inten- 
sity of the O2 lines, taken to be optically thin and produced in LTE, suggests that the 
temperature of the region is between 65 K and 120 K, and that the column density is 
N(02) = 6.5±lxl0^^ cm-^ Observations of J = 7 - 6 line obtained simultane- 
ously with those of O2 indicate emission in the same velocity range as observed for O2. 
Based on the C^^O emission at v > 10 km s~^ and assuming X(CO) = lxlO~^, we 
derive a H2 column density equal to 6.7x10^^ cm~^ and a fractional abundance X(02) 
~ 10"^ 

2. The hue centroid velocity is close to that of a number of species in a condensation ~ 6" 
west of IRc2, called Peak A, the Western Clump, or MF4 by various authors. Molecules 
including HDO, HC3N, NO, and HCOOCH3 have peak emission at velocities between 
10 km s"^ and 11.4 km s^^ towards this condensation, very similar to the observed 
velocities of the O2 lines. This is also the Cnt D localized maximum in the dust 
continuum emission. 

3. With only one pointing direction, we do not have definitive information on the exact 
source of the emission. The position of Peak A is ^ 21" from the direction of the 
present observations. The differential pointing of the two HlFl polarizations suggests 
that Peak A is more likely the source of the emission than the Hot Core, which is 
farther away from the beam pointing direction. 

4. If the 5" to 10"- diameter Peak A is the source of the observed O2 emission, the relative 
line intensities require that the source be quite warm, with T > 180 K, and the source 
column density N(02) = 4.6x10^^ cm~^. If the aforementioned C^^O emission also 
originates in this clump, our C^^O data imply an H2 column density equal to 3.7x10^^ 
cm~^. This is very similar to the H2 column density derived from interferometric 
observations of the dust in this clump, N(H2) = 5x10^^ cm"^. The resulting fractional 
abundance ix X(02) = 1.3x10"^, which is readily obtained in pure gas-phase chemistry. 

5. Explanations for the relatively large O2 fractional abundance compared to the results 
found by SWAS and Odin sateUites include (i) increase in gas-phase O2 as a result 
of warm grains having desorbed water ice mantles with subsequent reestablishment of 
standard gas-phase chemistry and (ii) enhancement of X(02) by shocks. 

(i) Heating dust grains to the ~ 100 K temperature required to desorb water mantles 
may have resulted from radiative heating, plausibly by nearby luminous sources. In this 
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case, the observed O2 lines could be emission from Peak A (possibly with contributions 
from other condensations that are all close enough to IRc2, IRc7, and other heating 
sources to have T grain ^ 100 K). The fractional abundance of O2 for this model implied 
by our observations is 1.3x10^^, which is more than an order of magnitude below that 
achievable. The time required to achieve the observed fractional abundance of O2 is 
moderate, < few x 10^ yr even for low cosmic ray ionization rate and varying inversely 
with this rate. Since desorption of water mantles is rapid, the total time for even the 
first high-mass star to produce a fractional abundance consistent with that observed, 
is compatible with the lifetime of high-mass stars. 

(u) A 10 km s~^ to 12 km shock can produce the observed O2 column density 
only if the emission is not significantly beam diluted. If the shock propagation di- 
rection is largely in the plane of the sky, the velocity centroid and line width may 
be consistent with what is observed. Shocks are generally assumed to be responsible 
for the widespread H2 emission observed in this region, but generally taken to have 
much higher velocities. J-shocks may produce UV radiation that maintains a signifi- 
cant abundance of atomic oxygen. This is required to be present in the gas which is 
then subjected to lower- velocity shocks which are optimum for producing the observed 
O2 column density. Shocks may also explain the water maser emission observed with 
Herschel. In this scenario, the O2 emission would have to come from numerous clumps 
distributed throughout the Herschel beams, possibly the same condensations respon- 
sible for the highly inhomogeneous H2 emission revealed by high angular resolution 
observations. 

The O2 velocities and widths, association with the high-column density warm Peak 
A source, and the plausibility of the warm dust model suggest that this compact 
condensation is responsible for the O2 emission we have observed. 

6. The resolution of the origin of enhanced O2 abundance will require both additional 
modeling and further observations. Our data do not suggest that the O2 abundance is 
sufficiently large to resolve the "oxygen problem" , referring to the question of allocat- 
ing the available oxygen among different forms (grains, ices, O, H2O, CO, and other 
molecules; Whittct 2010). Future observations will hopefully yield a more complete 
understanding of the astro chemistry of O2 and the distribution of oxygen in dense 
molecular clouds. 
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Appendix 1. HIFI Radiometer Noise Performance 



These relatively long integrations confirm the excellent performance of the HIFI instru- 
ment. The noise in the spectra of all three O2 spectra is in general quite close to that expected 
from the radiometer equation. The major limitation on the precision with which this can 
be determined is the possible presence of very weak spectral lines and of variable baseline 
ripple. We have carried out an extensive analysis of the 774 GHz data, for which the O2 line 
is accompanied by a moderately large range of frequencies without expected significant lines. 
We have data taken in March 2010 and September 2010 with a total of 26,848 s integration 
time. The noise (or fluctuation) bandwidth of the channels of the wideband spectrometer 
(WBS) used to obtain these data is determined by the channel frequency response, which is a 
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quasi-Lorentzian function (Tolls et al. 2004). The resulting noise bandwidth is Sunoise — 1-62 
MHz and the nominal channel spacing is 0.57 MHz. This will vary slightly from the intrinsic 
spectrometer value due to shifting of channels in the frequency calibration procedure. 

To clarify the comparison between theory and measurements, we have boxcar smoothed 
the data to a bandwidth of 2 MHz. According to the simple formula for the noise bandwidth 
for combination of channels (Ossenkopf 2008; Schieder & Kramer 2001) the noise bandwidth 
in this case should be close to 3 MHz. We utilize the single-sideband extracted spectra as 
described above, giving us a theoretical rms noise for a switched observation aiTMs) = 
ATnoise ( D SB ) ejj^ I \J Si'noisetiri t ■ whcic T{DSB) is thc double sideband noise temperature, 
5vnoise IS the uoisc bandwidth, and tint is the total (on + off) integration time. We use for 
all the data emb = 0.75, T{DSB) = 179 K for the V polarization channel, and T{DSB) = 
192 K for the H polarization channel. The noise temperatures varied by only a few K over 
the six-month interval spanned by the observations. 

We obtain for the V data a — 0.0062 K, and for the H data a — 0.0053 K, compared to 
theoretical values 0.0034 K and 0.0036 K, respectively. The ratio of observed to theoretical 
noise is 1.82 for V and 1.47 for H, with an average value of 1.65. These results are impacted 
by the difficulty of accurately determining the rms from data due to the limited number 
of "clean" channels and time dependent baselines, but they do indicate that the system 
performance is reasonably close to theoretical. 

The correctness of the noise bandwidth is confirmed by comparison of data smoothed 
to 4 MHz channel width (5 MHz noise bandwidth) with the 2 MHz channel bandwidth (3 
MHz noise bandwidth) data; the ratios for the different data sets have a mean value of 
0.76 compared to 0.78 (^3/5) expected. The rms of the 4 MHz smoothed data is 0.0050 
K for thc V channel and 0.0036 K for the H channel, which average to a factor 1.61 times 
greater than thc theoretical value. These results suggest that combining spectral channels 
reduces the noise in a manner close to that theoretically expected, and that the overall level 
of fluctuations is a factor ~ 1.6 greater than that expected strictly from the radiometer 
equation for integration times up to 2.7x10^ s. 

Appendix 2. Effect of Beam Dilution on Line Intensity Ratios 

In general, the integrated antenna temperature, /, can be written 



where in addition to symbols defined following equation 1, we have the antenna effective 
area A^i the normalized antenna response pattern P„, and the normalized column density 




(A2-1) 
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distribution, C„. For an extended uniform (C„ = 1) source filling the entire antenna pattern, 
the integral is the antenna solid angle, Qa- If the source is uniform but extended only to fill 
the main beam, the integral is the main beam solid angle, flMB- The main beam efficiency 
€mb is defined as emb — ^mb/^a- 

We can still use equation A2-1 by considering it to refer to a source that covers at 
most only the main beam, and a main beam-corrected integrated antenna temperature, 
Imb = / TAdv. From the antenna theorem, AgflA — A^, for a given antenna. Then we 
can rewrite equation A2-1 as 



Imb = / P„(Q)C„(Q)dQ , (A2-2) 

»7rfc \ImB J MB 

where the integral extends over the main beam. For a uniform source filling the main beam, 
the integral is just fluB and the integrated main beam temperature, I^^ is equal to the 
upper level column density multiplied by some constants. 



_AaNJic^ 

^ " 8nk ' ^^^""^^ 
where we use the superscript src to indicate that this intensity is essentially a characteristic 
of the source itself rather than of the antenna being used for the observations. For a uniform 
small source located in the direction of the antenna boresight , C„ = 1 within the extent 
of the source and = outside. The integral, which is nonzero over a range that is much 
smaller than the main beam size so P„ = 1, is equal to the source solid angle, Qs- Thus, the 
observed intensity, denoted by the superscript obs is given by 



jobs ^ AuiNuhcX^ Qs ,^2 4) 
STTk VLmb ' 

The ratio of the observed intensity to that of the source is given by 

R = 7°^7r- = ns/^MB . (A2-5) 

Using the antenna theorem again, and defining the aperture efficiency in terms of the physical 
area of the antenna through Ag = SA^p we find 

R = -^ApQsX-^ . (A2-6) 

^MB 



A low frequency (larger wavelength) transition has a small value of R, which can be thought 
of as arising from the small filling factor at the low frequency where the beam size is larger. 
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To analyze our Herschel observations of O2 in a reasonably simple manner, we consider a 
"small" source to be one that is smaller than the beam sizes at all three observing frequencies. 
Assuming a constant ratio of aperture to main beam efficiency, which is quite good for 
Herschel (Olberg 2010), R for two different transitions, 2 and 1, can be written as 

R{2/1) = [A(l)/A(2)]^ . (A2-7) 
Copyright 2011. All rights reserved. 
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